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The Aeroelastic Prediction Workshop brought together an international community of 
computational fluid dynamicists as a step in defining the state of the art in computational 
aeroelasticity. This workshop’s technical focus was prediction of unsteady pressure distri- 
butions resulting from forced motion, benchmarking the results first using unforced system 
data. The most challenging aspects of the physics were identified as capturing oscillatory 
shock behavior, dynamic shock-induced separated flow and tunnel wall boundary layer in- 
fluences. The majority of the participants used unsteady Reynolds-averaged Navier Stokes 
codes. These codes were exercised at transonic Mach numbers for three configurations and 
comparisons were made with existing experimental data. Substantial variations were ob- 
served among the computational solutions as well as differences relative to the experimental 
data. Contributing issues to these differences include wall effects and wall modeling, non- 
standardized convergence criteria, inclusion of static aeroelastic deflection, methodology 
for oscillatory solutions, post-processing methods. Contributing issues pertaining princi- 
pally to the experimental data sets include the position of the model relative to the tunnel 
wall, splitter plate size, wind tunnel expansion slot configuration, spacing and location of 
pressure instrumentation, and data processing methods. 
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I. Introduction 


UT t’s more an art than a science.” As engineers, we spend our careers creating methods, approaches, 
_L software and hardware to understand, explain and modify the physics that dominate a given problem. 
Somewhere between the imagining of a concept and the commercialization of a well-understood product is 
a long process that turns art into science. The AIAA Aeroelastic Prediction Workshop (AePW) and similar 
assessment activities are attempts to check in on where we stand as engineers in pursuing that process for 
different, and admittedly narrow, aspects of modeling aerospace vehicles. 

The AePW was held April 21-22, 2012 in Honolulu Hawaii, in association with the 53 rd AIAA/ASME 
/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference, and sponsored by the AIAA 
Structural Dynamics Technical Committee (SDTC). The workshop website contains the presentations from 
this meeting and other information required to begin analyses of the workshop configurations. A data 
summary of the AePW data sets has been presented previously. 2 

The AePW was a collaborative effort, in which time-accurate aerodynamic and aeroelastic computations 
using a vast array of codes and grids were exercised. The workshop analysts performed computations 
for their choices from among three configurations: the Rectangular Supercritical Wing (RSW) model; the 
Benchmark Supercritical Wing (BSCW) model; and the High REynolds Number AeroStructural Dynamics 
(HIRENASD) model. For each configuration, unforced system (steady) and forced oscillation (unsteady) 
analyses were performed and compared with existing experimental data sets. The unforced systems for the 
different configurations were treated either as rigid or flexible and analyzed with corresponding methods, as 
shown in table 1. The forced systems were similarly treated, analyzed with time-accurate solutions for the 
unsteady flow fields. Assumed-rigid and aeroelastic systems were analyzed differently. 

The RSW has simple geometry, a rectangular planform wing with transition strips on both the upper 
and lower surfaces near the wing leading edge. The RSW was tested at a fixed angle of attack and under 
forced oscillation pitching motion at frequencies well-separated from the structural dynamic modes. It was 
thus treated as a rigid model for the AePW. The test conditions examined for the AePW produced a flow 
field that was characterized as attached fully turbulent flow, with a moderate-strength shock. 

The BSCW has a similar simple geometry, again with a rectangular planform, transition strips and an 
assumed rigid structure. Fixed angle of attack and forced pitching oscillation data were also acquired for 
this model. The test conditions examined for the AePW corresponded to those that were observed during 
the experiment to generate strong shocks and have mixed attached/separated flow. 

The HIRENASD is an aft-swept aeroelastic wing with a generic fuselage fairing. This model was oscillated 
near the structural dynamic modal frequencies using wing root actuators. Inclusion of the HIRENASD ex- 
tends the choice of configurations to include more airplane-like systems and also introduces static aeroelastic 
coupling and very weak modal coupling into the problem. 

Table 1. Solution processes for AePW configurations 
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The AePW and its follow-on activities are collaborations within the aeroelastic community to assess the 
state-of-the-art in computational aeroelasticity, assess the experimental data available for performing this 
assessment and provide a roadmap forward. The intention is that the roadmap will address the required 
developments in both computational and experimental methods. 

This paper is intended to address the background, logistics, and reference information for the AePW 
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configurations and submitted computational results. Data and results from the workshop will be presented 
in companion papers. The remainder of this paper is structured addressing the motivation, genesis 

and strategy for this activity, followed by details of the technical tasks in preparation and execution of the 
workshop. 


II. Motivation 

The objective in conducting this workshop on aeroelastic prediction is to assess state-of-the-art Compu- 
tational Aeroelasticity (CAE) methods as practical tools for the prediction of static and dynamic aeroelastic 
phenomena. No comprehensive aeroelastic benchmarking validation standard currently exists, greatly hin- 
dering validation and state-of-the-art assessment objectives. 

Participants in the Drag Prediction Workshop (DPW) series 8 ’ 9 ’ ! and the High Lift Prediction Workshop 
(HiLiftPW) 1 1 ’ have performed quantitative and qualitative assessments of a significant cross-section of 
computational methods. These efforts have been used to determine the level of confidence that can be 
placed in computational results, focusing on steady-state rigid configurations. The AePW has been crafted 
to follow in the footsteps of these prior workshops, extending assessment efforts to time-accurate unsteady 
computational methods and coupled fluid/structure methods. 

The subsonic flow range properties are thought to be well-predicted by current methods, but the definition 
of “well-predicted” and the Mach number range for “good enough” agreement are phenomenon-dependent 
as well as end- usage-dependent. This workshop series aspires to provide a forum for unbiased and quantita- 
tive evaluation of different types of complexities, including flow phenomena, modeling, aeroelastic coupling 
strength, and configuration. 


III. Genesis of the Workshop 

The origins of this workshop lie within the NASA Fundamental Aeronautics Program. In 2009, the 
Subsonic Fixed Wing (SFW) project requested that an assessment of the state of the art in CAE methods 
be provided, assuming that a summary paper could be written- shown occurring in April 2009 on the timeline 
in figure 1. Comparable requests had been extended to those working steady-state computations of rigid 
vehicles. The response from that group was a summary of the extensive Drag Prediction Workshop (DPW) 
series. The comparatively sparse information regarding the state of the art in computational aeroelasticity 
sparked those involved to develop a plan that could better address the question of where CAE methods 
stood. After a year of groundwork, it was recognized that there might be interest throughout the technical 
community in participating in this assessment. The NATO Research and Technology Office (RTO) Aerospace 
Vehicle Technology (AVT) panel had an exploratory team developing a plan to perform validation wind 
tunnel tests. The missions seemed compatible and the AVT panel task group was queried about potential 
interest in collaborating on a workshop assessment, investigating existing experimental data sets. The DLR 
representative on the task group presented information regarding the HIRENASD configuration at the April 
2010 AVT meeting, suggesting it as a promising candidate for the collaboration. 

Thus, a collaboration was formed. A subsequent working group meeting was held in September 2010 
between RWTH Aachen — the organization that led the HIRENASD testing and research — and members of 
the NASA Langley Aeroelasticity Branch, working the SFW request. In October 2010, the official proposal of 
the NATO RTO AVT exploratory team changed from a new wind tunnel test to a joint exercise in aeroelastic 
prediction, utilizing existing databases. This was viewed by the RTO team as the logical first step, required 
before the requirements of a validation test could be properly identified, and viewed as a good path towards 
bolstering the technical case for performing a new set of validation experiments. 

A workshop exploratory committee was formed and requested participation from the aeroelastic com- 
munity, principally through the industry-led Aerospace Flutter and Dynamics Council (AFDC) meetings. 
A workshop Organizing Committee (OC) of volunteers was formed in December 2010. Planning, organiza- 
tion and conduct of the workshop was accomplished by the diverse group of international aeroelasticity and 
computational experts serving as the AePW OC, acknowledged in table 2. 

The AePW kick-off meeting was held in conjunction with the 2011 International Forum on Aeroelasticity 
& Structural Dynamics (IFASD). The OC solicited support from the AIAA SDTC, which then interacted 
with the AIAA to initiate the logistical matters for conducting the workshop. 

The focus of the OC was split between the technical aspects of the workshop (choosing the configurations 
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Figure 1. Timeline of events related to AePW. 


Table 2. AePW Organizing Committee members 


Name 

Affiliation 

Nation 

Kumar Bhatia 

Boeing Commercial Aircraft 

USA 

Josef Ballmann 

Aachen University 

Germany 

Eric Blades 

ATA Engineering, Inc. 

USA 

Alexander Boucke 

Aachen University 

Germany 

Pawel Chwalowski 

NASA 

USA 

Guido Dietz 

European Transonic Windtunnel (ETW) 

Germany 

Earl Dowell 

Duke University 

USA 

Jennifer Florance 

NASA 

USA 

Thorsten Hansen 

ANSYS Germany GmbH 

Germany 

Jennifer Heeg 

NASA 

USA 

Mori Mani 

Boeing Research & Technology 

USA 

Dimitri Mavriplis 

University of Wyoming 

USA 

Boyd Perry III 

NASA 

USA 

Markus Ritter 

DLR 

Germany 

David Schuster 

NASA 

USA 

Marilyn Smith 

Georgia Institute of Technology 

USA 

Paul Taylor 

Gulfstream Aerospace 

USA 

Brent Whiting 

Boeing Research & Technology 

USA 

Carol Wieseman 

NASA 

USA 
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and test conditions and preparing the information) and the organizational aspects of the workshop (recruiting 
analysts and getting information to and from them). From December 2010 until the kick-off meeting in June 
2011, the OC investigated available configurations and data sets. Initial discussions centered on which 
configuration should be chosen to achieve a balance between being of sufficient technical challenge and being 
sufficiently simple to provide a definitive benchmark. The efforts of organizing committee members also 
include generating grids for the workshop analysts to use, conducting baseline analyses to assess suitability 
of models and data sets and analyzing the experimental data sets. 

The AePW was held April 20-21, 2012, the weekend prior to the AIAA SDM conference. The agenda for 
the meeting is shown in Appendix A. 

IV. Building Block Strategy 

The AePW was a benchmarking activity where we drew on the logic and process that guides validation 
efforts. The AIAA Committee on Standards for Computational Fluid Dynamics’ written guide for verification 
and validation provides useful perspective. 

The fundamental strategy of validation is the identification and quantification of error and 
uncertainty in conceptual and computational models. The recommended validation method is to 
employ a building-block approach. This approach divides the complex engineering system of inter- 
est into . . . progressively simpler phases . . . The strategy in this approach is the assessment of how 
accurately the computational results compare with experimental data , with quantified uncertainty 
estimates, at multiple levels of complexity. Each phase of the process represents a different level 
of flow physics coupling and geometrical complexity. 

The OC followed these recommendations, dividing the complex problem of nonlinear unsteady aeroelas- 
tic analysis of an aerospace vehicle into simpler components. The components, or building blocks, were 
formulated to focus on specific aspects of the underlying physics. 

The classical building blocks of aeroelasticity are: 1) unsteady aerodynamics; 2) structural dynamics; and 
3) coupling between the fluid and the structure. In this first workshop, the OC chose to focus primarily on 
validating unsteady aerodynamic models and methods, with an initial venture into weakly coupled aeroelastic 
models. 

Within the category of unsteady aerodynamics, the problem was further refined through the selecting of 
configurations and test conditions. Decisions were made based on several criteria. The first criterion applied 
was the existence of a compatible and sufficient experimental data set. The second criterion applied for the 
initial workshop effort was perceived simplicity — simplicity of both configuration and aerodynanic phenom- 
ena. Phenomenologically, the OC chose to begin with moderately simple flow fields and moderately simple 
geometries and structures. In retrospect, additional consideration should have been given to simplifying the 
matrix of potential numerical solution choices. 

In this workshop, the OC chose to focus on transonic conditions for several reasons. Transonic conditions 
are often considered to be the most critical conditions with regard to aeroelastic phenomena. ’ ’ In the 
transonic range, various flow phenomena can initiate and produce severe aeroelastic issues such as flutter, 
limit cycle oscillation or buffet. As such, the most significant disagreements among computational results and 
between experiments and computations are observed. Coupling the criticality and the historically observed 
discrepancies in the transonic range drew the organizing committee to consider transonic predictions as the 
necessary starting point for discussion of workshop configurations and cases. 

In making the choice to examine transonic flow simplicity was undermined. Benchmarking first against a 
benign subsonic test condition seems an obvious requirement in retrospect, but in formulating the workshop, 
it was not a popular choice. The OC walked a tightrope between making the test problem interesting to 
analysts and making the test problem simple. Within the realm of transonic flows, the organizing committee 
attempted to make simplifying choices. 

There is an extensive range of unsteady transonic flow physics that could be considered and broken 
into building blocks. These phenomena include shocks of varying strength and position, as well as sepa- 
rated flow regions. The test cases were chosen with the idea of including cases with relatively simple flow 
characteristics — weak shock and attached flow — progressing to test cases with strong shocks and separated 
flow. Ideally, this test case progression would provide benchmarking and confidence in the ability to predict 
phenomena of increasing complexity. 
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Cases were chosen such that the flow field could reasonably be considered to be fully turbulent. The 
approach taken was to select cases where boundary layer trip strips artificially forced transition — eliminating 
transition location as a variable in the experimental data — with the goal of making the comparison with 
fully turbulent computational solutions easier. 

The objective of the AePW was to assess the state of the art of computational aeroelasticity methods as 
practical tools for predicting static and dynamic aeroelastic phenomena. The OC chose, as a first step towards 
this goal, to focus on the predictive capabilities for calculating steady and forced oscillation aerodynamics 
in the transonic regime for cases with forced turbulent flow transition. 

V. Preparation 

The workshop organizers invited the computational aeroelasticity community to apply best practices 
and state-of-the-art methods to predict unsteady aerodynamic characteristics for rigid or weakly-coupled 
aeroelastic systems. Within this scope, the test cases were laid out in building blocks of increasing complexity. 
The test configurations and conditions were selected in an attempt to advance in complexity from fully 
turbulent with attached flow and weak shocks to transient separation conditions with strong shocks and 
significant interactions between these flow features. 

The workshop participants analyzed three configurations: RSW, BSCW and HIRENASD. For each 
configuration, unforced system (previously termed “steady” or “static” ) and forced system (previously termed 
“unsteady” or “dynamic”) analyses were performed and compared with experimental data sets. The new 
terms to designate these data sets are required because the results of the workshop show that the unforced 
system data contains considerable unsteady, time-dependent, dynamic components. 

The configurations, described below, differentiate themselves in terms of flow phenomena, geometric 
complexity, Reynolds number and test medium. The geometric reference quantities for each configuration 
are specified in table 3. 




Table 

3. Reference quantities 



RSW 

BSCW 

HIRENASD 

Reference chord 

cref 

24 in 

16 in 

0.3445 m 

Model span 

b 

48 in 

32 in 

1.28571 m 

Area 

A 

1152 in 2 

512 in 2 

0.3926 m 2 

Dynamic pressure 

q 

108.9 psf 

200 psf 

40.055 kPa (for low Reynolds number cases) 
88.697 kPa (for high Reynolds number case) 

Moment 

X 

11.04 in 

4.8 in 

0.252 m 

reference 

y 

0 

0 

-0.610 nr 

point 

Z 

0 

0 

0 

Transfer 




Vertical 

function 


Pitch 

Pitch 

displacement 

reference 


angle 

angle 

(at x=l. 24521 m, 

quantity 




y=0. 87303 m) 

Test medium 


R-12 

R-134a 

Nitrogen 


V.A. Rectangular Supercritical Wing, RSW 

The RSW was chosen as the first configuration for the workshop due to its simple geometry and structure, 
the relatively benign transonic flow characteristics and the availability of forced oscillation data. The RSW 
test cases, listed in table 4, were chosen to focus on the steady and unsteady aerodynamic solutions and 
their variation. Data at Mach 0.825 was the highest test Mach number where experimental data existed with 
transition strips on the model. At this condition, the flow includes a moderate strength shock. Two angles 
of attack were chosen from among the unforced system data. The lower angle, 2°, was chosen to match the 
available oscillatory data’s mean angle of attack. At this condition, a moderate strength shock was observed 
with some potential for shock-separated flow. 
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Table 4. RSW analysis conditions for AePW. 


Mach 

Mean a 
(deg) 

Pitch Osc. Freq. 
(f, Hz) 

Pitch Osc. Amp. 
(6>, deg) 

Reduced Freq. 
w c / 2Voo 

Re c 

no 6 

0.825 

2.0 

0 

0.0 

0 

4.0 

0.825 

4.0 

0 

0.0 

0 

4.0 

0.825 

2.0 

10 

1.0 

0.152 

4.0 

0.825 

2.0 

20 

1.0 

0.304 

4.0 


The simple geometry and structure of this wing were important factors in the decision to include it as the 
first test configuration. It has a rectangular planform wing as shown in figure 2, with a 12% thick supercritical 
airfoil. Adding to the simplicity, the model was designed with the goal of being structurally rigid so that 
there was negligible static aeroelastic deflection under the dynamic pressure loads applied in the testing and 
the aeroelastic modes had natural frequencies sufficiently above the frequency of applied oscillations. The 
first bending mode natural freqency of the RSW was measured at 34.8 Hz. Chordwise rows of unsteady 
pressure sensors were installed at 4 span stations, providing the principal data for this configuration. The 
model was tested in the NASA Langley Transonic Dynamics Tunnel (TDT) in 1982, 1 ' with R-12 freon as 
the test medium. It is shown in its test configuration in figure 3, mounted on a small splitter plate offset 6 
inches from the wind tunnel wall. The wall slots of the were open for the data acquired during RSW testing. 
The mounting geometry and the tunnel configuration play important roles in complications that arose in 
analyzing the RSW. These testing issues are explored in more depth in a companion paper. ' 




Figure 3. RSW, mounted in the wind tunnel. 


Figure 2. RSW model geometry and layout. 


Time-averaged (mean) data was obtained from data sets where the model held at a fixed angle of attack. 
These mean data sets were used to characterize the expected steady behavior of the configuration. Dynamic 
data was acquired by oscillating the model in a pitching motion about the 46% chord line. Two unforced 
system cases, 2° and 4° angle of attack, allowed evaluation of the angle of attack effects from a steady- 
state perspective. The two forced pitch oscillation frequencies (10 and 20 Hz) allowed evaluation of the 
methods abilities to distinguish frequency effects. Frequency separation between the structural dynamic 
modes and the excitation forcing frequency was used to minimize the coupling of the aerodynamics and the 
structural dynamics. The non-zero mean angle of attack introduces a wing loading bias for which code-to- 
code comparisons can be made. 

The steady and unsteady data was originally published as a NASA Technical Memorandum and later 
included in an appendix to AGARD and RTO reports. The RSW experimental data has been widely 
published. Unfortunately, the only available data is the published data. It consists of mean pressures for static 
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cases and complex pressures for forced oscillation runs. The model is unavailable for retesting or inspection, 
although the splitter plate was located and its geometry was measured and used in grid generation. 

The principal deficiency of the RSW configuration is the proximity of the wing relative to the wind tunnel 
wall, which effectively engulfs the inboard portion in the wall boundary layer. This issue and the implications 
are thoroughly discussed in a companion paper. 

V.B. Benchmark Supercritical Wing, BSCW 

The BSCW was chosen as the second workshop configuration due to its simple geometry and structure, 
more complex transonic flow characteristics and the availability of forced oscillation data. The BSCW has a 
rectangular planform as shown in figure 4, with a supercritical airfoil, SC(2)-04114. The BSCW was treated 
as a rigid wing in the AePW computations. The wing was tested at frequencies that encompass the first 
bending mode frequency, but only data at excitation frequencies well below this value were used for the 
AePW. The spanwise first bending mode has a frequency of 24.1 Hz, the in-plane first bending mode has 
a frequency of 27.0 Hz and the first torsion mode has a frequency of 79.9 Hz. Any static aeroelastic effect 
was assumed to affect the wing uniformly, and was accounted for by using the measured wing root angle of 
attack in the analysis conditions. The BSCW, shown in figure 5, was tested in the TDT in an R-134a test 
medium. The model was mounted on a large splitter plate, sufficiently offset from the wind tunnel wall (40 
inches) to be outside of the tunnel wall boundary layer. The testing was conducted with the sidewall slots 
open. The models instrumentation is limited to one row of 40 in-situ unsteady pressure transducers at the 
60% span station. Boundary layer transition was fixed at 7.5% chord using size 30 grit. All data to be used 
in this study was obtained at Mach 0.85 and a dynamic pressure of 200 psf, fixing the Reynolds number at 
4.49 million based on wing chord. The AePW BSCW test cases are summarized in table 5. 


T 


16 " 


■y 


Pitch axis, 
30% chord 


32 " 


Transition strip, 
7.5% chord 


' f x 


Pressure 
transducers, 
60% chord 


Figure 4. BSCW model geometry and layout. 


Figure 5. BSCW, mounted in the wind tunnel. 



Table 5. BSCW analysis conditions for AePW 


Mach 

Mean a 
(deg) 

Pitch Osc. Freq. 
(f, Hz) 

Pitch Osc. Amp. 
{0, deg) 

Reduced Freq. 
cue ! 2 Vqo 

Re c 

no 6 

0.85 

5.0 

0 

0 

0 

4.49 

0.85 

5.0 

1 

1.0 

0.009 

4.49 

0.85 

5.0 

10 

1.0 

0.090 

4.49 


Dynamic data was obtained for the BSCW by oscillating the model in a pitching motion about the 30% 
chord using an oscillating turntable. Unforced system quantities for this configuration were calculated as 
the time-averaged (mean) values obtained without excitation. To better capture the information, maximum 
and minimum values that occurred during the unforced system time records were also included in the 
comparison data set. 
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From the standpoint of aerodynamic complexity, this configuration was chosen because the experiment 
exhibited nonlinear responses to the sinusoidal pitch excitations, pointing to shock-separated transient flow. 
While there are fewer pressure measurements than for the RSW configuration, the time history data records 
are available for all test conditions. The model is also available for inspection and retesting if desired. The 
engineers who conducted the testing and their detailed test notes are also available for consultation regarding 
the model, test conditions, data content and other insights. 

V.C. High REynolds Number AeroStructural Dynamics, HIRENASD 

The HIRENASD was chosen as the final configuration for the workshop because of it’s more airplane- 
like geometry, weak aeroelastic coupling and thorough transonic data set including forced modal-frequency 
oscillations.' ’ u ’ ’ The HIRENASD was tested in the European Transonic Wind tunnel (ETW) in 
2007, shown in figure 6 using Nitrogen as the test medium. The model has a 34° aft-swept, tapered clean wing, 
with a BAG 3-11 supercritical airfoil profile. The test article is a semi-span model, ceiling-mounted through 
a non-contacting fuselage fairing to a turntable, balance and excitation system, shown in figure 7. The model 
and balance were designed to be very stiff, with well-separated modes. The first two bending modes have 
frequencies of approximately 27 and 79 Hz; the first torsion mode has a frequency of approximately 265 Hz. 
The HIRENASD forced motion test cases used for the AePW are oscillated at the second bending mode 
frequency using piezoelectric stacks located inside the balance/standoff structure at the wing root. The 
models instrumentation includes 259 in-situ unsteady pressure transducers at 7 span stations. In addition 
to the unsteady pressures, balance measurements and accelerations were obtained. For a small set of data 
points, wing displacements were also extracted via stereo pattern tracking. 

For the workshop comparisons, data from a single test point was used to provide both steady and unsteady 
data. Data sets with forced oscillations were obtained during the wind tunnel test by sending a low amplitude 
sinusoidal command, followed by a “rest period”, then followed by a high amplitude sinusoidal command. 
A subset of the data obtained during the rest period is used to calculate the unforced or steady comparison 
data. The forced oscillation data was obtained by differential forcing at a specified modal frequency. 

The HIRENASD test conditions used for the AePW are presented in table 6. Two Reynolds numbers, 7.0 
million and 23.5 million based on reference chord, were analyzed for Mach 0.8. The lower Reynolds number 
case has an angle of attack of 1.5°, while a more challenging angle of attack of -1.34° (corresponding to the 
zero-lift condition) has been selected for analysis at the higher Reynolds number. A third test case was added 
at Mach 0.7 to provide a benign case that could presumably be analyzed without substantial problems. For 
the low Reynolds number workshop test case conditions, boundary layer transition strips were affixed to the 
upper wing surface at 12-15% chord and to the lower surface at 5% chord; the transition strips were removed 
for testing at the 23.5 million Reynolds number condition. 


Table 6. HIRENASD analysis conditions for AePW 


Mach 

Mean a 
(deg) 

Exc. Freq. 
(f, Hz) 

Osc. Amp. 
(Az, mm) 

Reduced Freq. 

w c / 21/*, 

Re c 

*10 6 

0.70 

1.5 

0 

0 

0 

7 

0.80 

1.5 

0 

0 

0 

7 

0.80 

-1.34 

0 

0 

0 

23.5 

0.70 

1.5 

79.3 

2.0 

0.34 

7 

0.80 

1.5 

78.9 

2.4 

0.39 

7 

0.80 

-1.34 

80.4 

0.9 

0.38 

23.5 


The HIRENASD model was chosen as the initial coupled aeroelastic analysis configuration. The high 
stiffness and modal spacing produces weak aeroelastic coupling, making it a good entry-level basis of eval- 
uation. One of the challenges in using the HIRENASD was generating a structural dynamic model that 
was well-understood by the workshop team. Variations in the bending mode frequency in publications, and 
the use of various available structural models had to be understood and resolved. The final model that was 
seledted for use included a substantial portion of the model mounting system hardware, as shown in figure 8. 
The resulting second bending mode shape is shown in figure 9. 
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Figure 6. HIRENASD, mounted in ETW. 


damping flange noncontact 

(balance connection) labynnth sealing 



Figure 7. HIRENASD model geometry and layout. 


10 of 21 


American Institute of Aeronautics and Astronautics 




Figure 8. Final HIRENASD finite element model. 



Figure 9. HIRENASD 2nd Bending mode shape. 


V.D. Requested comparison data 

For each configuration to be analyzed, results from three studies were requested: a convergence study, steady- 
state analysis and time-accurate response due to forced oscillations. The convergence studies that were 
requested were both spatial and temporal, performed separately. A grid convergence study was requested 
using a coarse, medium and fine grid as defined in the gridding guidelines to be discussed below. The analysts 
were asked to submit both steady-state and forced oscillation results for each of the grid refinements. The 
requested comparison parameters were the integrated load coefficients (lift coefficient, drag coefficient and 
pitching moment coefficient) for the steady-state calculations, and the frequency response functions due to 
reference excitation for the unsteady calculations. 

The steady-state computational results were asked to include the pressure coefficient distribution as 
a function of normalized chord location and normalized span station for each of the pressure transducer 
locations on the respective experimental models. Many of the analysis teams, however, provided their 
results at all of the grid locations at the span stations of interest. Those analysts examining the HIRENASD 
configuration were also asked to provide the vertical linear displacements and twist angles at span stations 
corresponding to the displacement sensor locations. 

Results from oscillatory analyses were asked to include the frequency response functions (FRFs) of the 
pressure coefficients due to displacement. These FRFs were requested to be calculated as magnitude and 
phase values for a single frequency corresponding to the frequency of the excitation. Time histories of the 
pressure coefficients were also requested at a chosen span station for each of the configurations; the chord 
location was specified for one of the locations, and the other chord location was to be generated at the 
location corresponding to the peak of the FRF. 

Definition of the reference displacement signal was different for the different configurations. The RSW 
and the BSCW are both oscillated in pitch; analyses of both configurations should use the angle of attack 
displacement as the excitation source in computing the frequency response functions. The reference signal 
for HIRENASD should be the vertical displacement at x=l. 24521 m, y=0. 873034 m, corresponding to the 
location of the wingtip accelerometer. 
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V.E. Geometry and grid preparations 

Model geometry files were made available for each configuration. Grids generated by members of the AePW 
OC were provided to analysts for optional use. The majority of the analysts, however, chose to generate 
their own grids. Gridding guidelines were established and are detailed in Heeg et al. These gridding 
guidelines for the Aeroelastic Prediction Workshop are adopted from the guidelines developed for the DPW 
and HiLiftPW. Mavriplis et al." summarize the grid-related lessons gleaned from the experiences of DPW. 
These guidelines remained relatively unchanged over the course of DPW II-IV and were maintained in 
the HiLiftPW; the guidelines codify much of the collective experience of the applied CFD community in 
aerodynamic grid generation practices for steady analyses. In DPW V, they have reinvigorated their studies 
of grid quality effects. The past workshop gridding guideline development efforts have not addressed issues 
particular to performing time-accurate calculation for unsteady problems. Thus, grid development guidelines 
directly pertinent for the AePW have not been established. For this first workshop in aeroelastic prediction, 
the steady grid guidelines used in DPW II-IV and HiLiftPW were adopted. 

The gridding guidelines included specification of the cell spacing normal to the walls, boundary layer cell 
growth rate, wing edge cell spacings, number of cells on the trailing edge face cells and fuselage cell size. 
The progression from coarse to fine grids required that the grid size increase by a factor of 3 for each level 
of refinement for unstructured grids, and by a factor of 1.5 for structured grids. Special effort is required 
to ensure that sequences of coarse, medium and fine meshes constitute a consistent family of grids suitable 
for a grid convergence study. This entails the preservation of mesh topology, stretching factors, and local 
variations in resolution as much as possible between grids of the same sequence. 

V.F. Database structuring and software development for information tracking and compar- 
isons 

Data was submitted by the analysts using spreadsheets with predefined templates, with separate templates 
for steady and unsteady cases. An macro was written to save the sheets as text files, with one sheet for 
each submitted case result. The information includes hie name mapping for tracking of data submittals and 
updates of information. A program was coded to read all of the existing text hies and generate data hies 
suitable for data processing software, named based on hie contents. Another program was written to combine 
the data hies for each conhguration and test case into a single data base. Numerous checks on descriptive 
helds were necessary to create data set identihers within the database to track text-based analyst-provided 
parameters. 

Data requested from the analysts included steady pressure coefficient distributions, including locations 
of those calculated pressures, magnitude and phase of frequency response functions at dominant response 
frequency for each case, static aeroelastic deflections for the HIRENASD case, and a few example time 
histories of pressure coefficients. The data was requested to correspond to the wind tunnel model sensor 
locations on the upper and lower surfaces of the wings. 

V.G. Examination of data: Comparison software 

Software was written to produce comparison plots of the data sets. Modular software was developed for quick 
database and plot updates as new datasets were provided. These modular codes provided the capability 
of sorting the data, treating any of the information helds in the database as the independent variables (for 
example turbulence model or grid refinement level) . Separate modules generated comparison plots of pressure 
coefficient as a function of chord location, and magnitude and phase of frequency response functions versus 
chord location. Integrated lift, drag and pitching moment coefficient plots are generated as functions of grid 
factor or time step, providing consolidated convergence information. For the HIRENASD cases, vertical 
dehections at span stations were also compared. 

VI. Execution 

Seventeen analysis teams from 10 nations participated by providing computational results for the work- 
shop, hgure 10. The analysis teams, organizations and conhgurations analyzed by each team are listed in 
table 7. The choice of how solvers and associated turbulence model and flux construction selections are listed 
by conhgurations in tables 8 (RSW), 9 (BSCW), and lO(HIRENASD). 
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The majority of the analysis teams utilized Reynolds Averaged Navier-Stokes (RANS) flow solvers. RANS 
equations are time-averaged equations of motion for fluid flow. The time-averaged assumptions involve sepa- 
ration of chaotic velocity Auctions from the mean Aow velocity. This in turn requires the use of a turbulence 
model which is usually tuned for a speciRc Aow physics situation. The unforced system analyses were in 
general performed by converging the RANS solutions to steady state. The forced oscillation simulations were 
performed using unsteady RANS (URANS) codes, solved in a time-accurate manner with subiterations to 
converge the solution at each of the time steps. 

Two HIRENASD analysis teams did not use RANS Aow solvers. As noted in table 10, HIRENASD 
analysis team I performed Euler Aow solutions using ZEUS code. HIRENASD analysis team L performed 
full-potential Aow solutions using the ST Aow solver code. 

It should also be noted that HIRENASD analysis team C performed direct-coupled simulations rather 
than implementing a structural modal solver. The modal solution approach was used by all other teams in 
analysis of the HIRENASD case. 

VI. A. Grids 

The workshop analysts overwhelmingly chose to 
generate their own grids. Of the 17 analysis teams, 
only 4 teams chose to use a grid generated by some- 
one else. The diversity of grid types and sizes uti- 
lized for each of the conRgurations are summarized 
in tables 11, 12 and 13. The number of grid points 
that deAned Coarse(C), Medium(M) and Fine(F) 
grids for each analysis team are given in these same 
tables and plotted by conRguration in Agures 11, 

12 and 13. The qualitative descriptors are those 
assigned by the analysis teams, rather than being 
based on the submitted grid size quantihers. It is 
noteworthy that some analysis teams coarse grids 
were as large as other analysis teams medium grids, 
particularly for the HIRENASD case. Also, BSCW analysis team C submitted results for two medium grids, 
creating the job shown in in Agure 12. 

The descriptions of each of these grids provides a basis for comparing the grids, however, the numbers do 
not represent the entire story. Among these grids are diRerences in regions of different element size gradients 
and outer boundary computational domain shape. 



Figure 10. Analysis teams from 10 nations participated 
in the first Aeroelastic Prediction Workshop; Experimen- 
tal data provided by RWTH Aachen, Germany &: NASA, 
USA. 



Figure 11. RSW grid resolutions. 
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Table 7. AePW analysts, affiliations &: cases analyzed 


Affiliation Analysis Team Members 

RSW 

BSCW 

HIRENASD 

NASA Pawel Chwalowski 

X 

X 

X 

ANSYS Germany GMBH Thorsten Hansen, 

Angela Lestari 

X 

X 

X 

University of Wyoming Dimitri Mavriplis, 

Mike Long, 

Zhi Yang, 

Jay Sitaraman 

X 

X 

X 

RUAG Aviation Alain Gefiri, 

Daniel Steiling 

X 

X 

X 

NASA David Schuster, 

Andrew Prosser 

X 

X 


Swedish Defense Research Agency, FOI Mats Dalenbring, 

Adam Jirasek 


X 

X 

Technion University IIT Daniella Raveh 


X 

X 

Georgia Institute of Technology Marilyn Smith, 

Benn Mann 


X 


University of Liverpool Sebastian Timme 

X 



NLR Bimo Pranata, 

Bart Eussen 
Jaap van Muijden 



X 

ONERA Anne-Sophie Sens, 

Jean-Pierre Grisval 



X 

DLR Markus Ritter 



X 

Istanbul Technical University Melike Nikbay, 

& Pinar Acar, 

Cagri Kilic, 

ZONA Technology, Inc. Zhichao Zhang 



X 

Politecnico de Milano Sergio Ricci, 

Andrea Parrinello, 
Giulio Romanclli 



X 

MSC and Jack Castro, 

Metacomp Beerinder Singh 



X 

Boeing Research & Technology Mori Mani, 

Andrew Cary, 
Larry Brase 



* 

X 

CD-Adapco Alain Mueller, 

Sergey Zhelzov 



* 

X 


Analyses performed were different from those required for comparison with other AePW datasets. 
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Table 8. AePW RSW flow solution information 


Analysis 

Team 

Software 

Name 

Turbulence 

Model* 

Flux 

Construction 

Flux 

Limiter 

Oscillatory 

Solution 

Method 

A 

NSMB 

SA 

Unknown 

None 

Elastic+TFI 

B 

FUN3D 

SA 

Roe 

Venkat 

Elastic 

C 

CFL3D 

SA 

Roe 

None 

Modal+TFI 

D 

ANSYS CFX 

SST 

2nd Order Upwind/ 
Rhie Chow 

Bartli & 
Jesperson 

Diffusion equation 

E 

NSU3D 

SA 

Matrix Artificial 
Dissipation 

None 

Full grid motion 

F 

PMBvl.5 

SA 

Oslier 

MUSCL+ 
van Albada 

Full grid motion 


Spalart-Allmaras (SA), Shear Stress Transport (SST) 


Table 9. AePW BSCW flow solver information 


Analysis 

Team 

Software 

Name 

Turbulence 

Model* 

Flux 

Construction 

Flux 

Limiter 

A 

NSMB 

SA 

Unknown 

None 

B 

FUN3D 

SA 

Roe 

Venkat 

C 

CFL3D 

SA 

Roe Flux 
difference splitting 

None 

D 

NSU3D 

SA 

Central difference with 
matrix dissipation 

Unknown 

E 

ANSYS CFX 

SST 

2 nd Order Upwind/ 
Rhie Chow 

Barth & Jesperson 


Spalart-Allmaras (SA), Shear Stress Transport (SST) 



Figure 12. BSCW grid resolutions. 
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Table 10. AePW HIRENASD flow solver information 


Analysis 

Team 

Software 

Name 

Turbulence 

Model* 

Flux 

Construction 

Flux 

Limiter 

A 

ENFLOW 

k-w 

Central difference with 
artificial dissipation 

TVD 

B 

NSMB 

k-w 

unknown 

None 

C 

CFD++\ NASTRAN 

k-e 

HLLC 

Compressive-MinMod 

D 

EZNSS 

SA 

HLLC 

Venkat 

E 

EDGE 

SA 

Central difference 1 

None 


EDGE 

SA 

Roe 2 

Venkat 

F 

TAU 

SA 

Central scheme 
with scalar 
dissipation 

None 

G 

elsA 

SA 

Jameson 

None 

H 

NSU3D 

SA 

Central difference with 
matrix dissipation 

None 

I 

ZEUS 1 ' 

G 

Central difference with 
JST artificial dissipation 

None 

J 

FUN3D 

SA 

Roe 

Venkat 

K 

ANSYS CFX 

SST 

2nd Order upwind/ 
Rhie Chow 

Barth & Jesperson 

L 

ST 1 

n/a 

unknown 

None 

M 

AeroFoam 

SA 

Roe + LW 

vanLeer 


Spalart-Allmaras (SA), Shear Stress Transport (SST), Green’s integral boundary layer method (G), variants of k-o; 
(k-cj), 2-equation realizable k-e (k-£) 
t Euler flow solver 
$ Full potential equation flow solver 

1 Denoted Edge flux construction 1 in plot legends 

2 Denoted Edge flux construction 2 in plot legends 


Table 11. AePW RSW submitted grids 


Analysis 

Grid 

Element 

Solver 

Number of Nodes or Cells, (millions) 

Team 

Type* 

Type 1 

Type 1 

Coarse 

Medium 

Fine 

A 

Str 

Hex 

Cell 

3.38 

9.91 

27.0 

B 

Unstr 

Mix 

Node 

2.88 

7.07 

18.23 

C 

Str 

Hex 

Cell 

0.18 

1.42 

11.18 

D 

Str 

Hex 

Node 

1.91 

5.89 

15.42 

E 

Unstr 

Mix 

Node 

2.87 

7.07 

18.28 

F 

SMB 

Hex 

Cell 

2.32 

6.60 

18.63 


* Structured (Str), Unstructured (Unstr), Structured MultiBlock (SMB) 
T Hexagonal (Hex), Mixed Hexagonal & Tetrahedral (Mix) 

$ Cell-centered (Cell), Node-centered (Node) 
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Table 12. AePW BSCW submitted grids 


Analysis 

Grid 

Element 

Solver 

Number of Nodes or 

Cells, (millions) 

Team 

Type* 

Type* 

Type* 

Coarse 

Medium 


Fine 

A 

Str 

Hex 

Cell 

3.79 

9.48 


30.32 

B 

Unstr 

Mix 

Node 

2.97 

9.01 


26.79 

C 

Str 

Hex 

Cell 

0.14 

1.07 


8.40 






1.59 



D 

Unstr 

Mix 

Node 

2.97 

9.01 



E 

Str 

Hex 

Node 

1.49 

5.03 


13.93 


* Structured (Str), Unstructured (Unstr), Structured MultiBlock (SMB) 
' Hexagonal (Hex), Mixed Hexagonal & Tetrahedral (Mix) 

$ Cell-centered (Cell), Node-centered (Node) 


Table 13. AePW HIRENASD submitted grids 


Analysis 

Grid 

Element 

Solver 

Number of Nodes or 

Cells, (millions) 

Team 

Type* 

Type* 

Type* 

Coarse 

Medium 


Fine 

A 

SMB 

Hex 

Cell 


10.66 



B 

Str 

Hex 

Cell 


9.69 



C 

Unstr 

Mix 

Cell 

2.93 

8.36 


24.18 

D 

SMB 

Hex 

Cell 


12.61 



E 

Unstr 

Mix 

Node 

6.50 





Unstr 

Mix 

Node 

6.36 




F 

Unstr 

Mix 

Node 

1.03 

2.45 


7.21 

G 

SMB 

Hex 

Cell 

7.20 




H 

Unstr 

Mix 

Node 

6.36 

19.06 



I 

Str 

Quad 

Cell 

0.56 




J 

Unstr 

Mix 

Node 

6.36 

19.06 


56.31 

K 

Str 

Hex 

Node 


10.03 



L 

Unstr 

Tet 

Node 

0.14 




M 

Unstr 

Mix 

Cell 

1.63 





* Structured (Str), Unstructured (Unstr), Structured MultiBlock (SMB) 
t Hexagonal (Hex), Mixed Hexagonal & Tetrahedral (Mix), Tetrahedral (Tet), Quadrilateral 
(Quad) 

$ Cell-centered (Cell), Node-centered (Node) 
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Figure 13. HIRENASD grid resolutions. 


VII. Concluding Remarks 

The AePW was envisioned as a step in assessing the state of the art in computational aeroelasticity. It 
has demonstrated that the current state of the art in flow solvers is RANS, with the perceived best practice 
turbulence model being Spalart-Allmaras. We make this statement based on choices of the analysts, rather 
than on any defining characteristic in the results. The AePW analysis teams generally had these types of 
codes available and sufficient expertise and resources to exercise them, at least for the geometric complexity 
being captured in the AePW configurations. 

This is in contrast with the perceived state-of-the-current-practice. The current practice for production 
aeroelastic analysis is coupled analysis of linear aerodynamic and linear structural dynamic models. The 
perception is that the linear analysis community did not provide computational results for the workshop 
because they understand the level of the errors that they are accepting with their calculations, and they 
understand the limitations that they are imposing on their results in terms of capturing the flow physics. 
There are correction methods for these linear results that range from public and simple to proprietary and 
highly configuration dependent. 

A state-of-the-art assessment, however, necessarily includes more than a description of the paintbrushes. 
Companion papers present analysis and discussion of the computational and experimental data generated 
by the AePW community. >4 , <>. 

In this paper, we have attempted to summarize the AePW from an organization standpoint, hoping that 
others will see how and why decisions were made and where those decisions directed this workshop. Figure 1 
shows the AePW follow-up events, including special sessions at technical venues. These events are being 
conducted to more formally publish the workshop analyses and findings. The first of these special sessions 
is the AIAA Aerospace Sciences Meeting, where the current paper is being presented. Additional special 
sessions will be conducted at the AIAA Structures, Structural Dynamics and Materials Conference in April 
2013 and the 2013 International Forum on Aeroelasticity and Structural Dynamics. At each special session, 
an open discussion is planned with the objectives being twofold: 1) a frank discussion of the interpretation 
and significance of the results to date, and 2) provide a forum for continued discussion on the path forward. 

At each progressive session, an overview paper will capture the content of the prior discussion sessions. 
The intention in organizing these special sessions is for the workshop participants to continue to collaborate 
and push towards understanding the results beyond their own individual calculations. 
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all retired gracefully. 

The OC members have brought their own perspective and objectives to this workshop activity. The 
opinions of this group guided the technical decisions and steered the preparation efforts. The aeroelasticity 
community has been known for its technical throwdowns. We are particularly grateful that this group was 
vocal, yet civilized. 

“Imitation is the sincerest form of flattery.” ' 3 " We are hopeful that the organizers and participants of the 
HiLiftPW and DPW don’t mind that we have borrowed liberally from the organizational aspects, websites 
and technical perspectives of these teams. We gratefully acknowledge that they paved the way for us. 

The NASA portion of this work was funded by the Fundamental Aeronautics Program’s Fixed Wing and 
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IX. Appendix A: AePW agenda 

The agenda for the workshop is given in table 14. 

Presentations can be found online at https://c3.nasa.gov/dashlink/projects/47/. 
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Table 14. Agenda for the Aeroelastic Prediction Workshop 


Topic 

Presenter 

Organization 

Welcome & Workshop Overview 

Jennifer Heeg 

NASA 

Experimental Data Reduction Methods 

Jennifer Heeg 

NASA 

Overview of Rectangular Supercritical Wing Test Case 

Boyd Perry 

NASA 

RSW Analysis Presentations 

Pawel Chwalowski 

NASA 

Session Chair: 

Thorsten Hansen 

ANSYS Germany GMBH 

Alexander Boucke, 

Dimitri Mavriplis 

University of Wyoming 

RWTH Aachen 

David Schuster 

NASA 


Daniel Steiling 

RUAG Schweiz AG 


Sebastian Timme 

University of Liverpool 

RSW Comparisons & Discussion 



Overview of Benchmark Supercritical Wing Test Case 

Rob Scott 

NASA 

BSCW Analysis Presentations 

Pawel Chwalowski 

NASA 

Session Chair: 

Thorsten Hansen 

ANSYS Germany GMBH 

Brent Whiting, 

Dimitri Mavriplis 

University of Wyoming 

Boeing 

David Schuster 

NASA 


Daniel Steiling 

RUAG Schweiz AG 


Marilyn Smith 

Georgia Tech 

BSCW Data Comparisons 



Overview of HIRENASD Test Case 

Alexander Boucke 

RWTH Aachen 

Structural Dynamics Modeling for HIRENASD 

Carol Wieseman 

NASA 

HIRENASD Analysis Presentations 

Daniel Steiling 

RUAG 

Session Chair: 

Bart Eussen 

NLR 

Kumar Bhatia, 

Dimitri Mavriplis 

University of Wyoming 

Boeing 

Markus Ritter 

DLR 


Thorsten Hansen 

Ansys 


Mats Dalenbring 

FOI 


Pawel Chwalowski 
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